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We study superconductivity in multilayer copper oxides, in the frame of a realistic microscopic 
formulation. Solving the full temperature dependent BCS gap equations, we obtain a maximum in 
the transition temperature Tc for M=3 or 4 Cu02 layers in the unit cell for appropriate values of 
the interlayer tunneling (negative pair tunneling), and via the consideration of the doping imbalance 
between the inner and outer layers. This is the ubiquitous experimental result for Ca intercalated 
copper oxides, as opposed to other intercalating elements. Further, using a restricted set of param- 
eters, we obtain an exact fit of rc(M=l-4) for five different Ca intercalated homologuous copper 
oxide families. 

PACS numbers: 74.72.-h, 74.20.Fg, 74.62.-c, 74.78. Fk 



A long standing puzzle of the high transition temperature Tc cuprate superconductors is the dependence of T^ 
^ on the number M of Cu02 layers in the unit cell. Experiments show ubiquously that within any Ca intercalated 
O ^ homologuous cuprate family, i.e. for materials sharing the same charge reservoir block but having varying M (here 
^ \ Ca is between the layers), T^ is maximum for M=3 - a summary of experiments appears in e.g. [l[ - with the M=4 
■ exception of the Tl2Ba2CaM-iCuM02M-i-3-i-i5 family 0], following a dome-type dependence on M. That is at optimum 
' doping for any M. As mentioned in other spacing elements, such as Y, Ba and Sr yield lower for the multilayer 
' materials compared to the single layer material. Setting aside the central issue of superconductivity within a single 
i' \ Cu02 layer, the issue of what determines Tc in a multilayer system is a very important one. Inter alia, it points to 
'■^ ■ the role of out-of-plane mechanisms and to how Tc can be maximized for airy given cuprate family. 
C \ This issue has been treated in a number of theoretical papers, as in [H. l4llo|. with limited success. E.g. in ref. 
^ ■ [H a Bardeen-Cooper-Schrieffer (BCS) model (different from ours) was considered, and a saturating increase of Tc 
as a function of M=l-oo was obtained. In ref. [loj the dome with a maximum for M=3 was obtained for the 
superconducting order parameter (not directly for Tc itself) in the frame of a Landau-Ginzburg functional treatment, 
and by appealing to the putative ddw order parameter. 

The cuprates fall into the realm of strong correlations, thus calling for an Eliashberg-type treatment 11- l^l- Given 



> 

\f-^ _ the lack of a definitive theory for the single Cu02 layer superconducting state, our BCS treatment is a decent attempt 
\^ • towards the understanding of the multilayer cuprates. We hereby consider a realistic microscopic model in the frame 
of the BCS formulation |14|, which does yield the experimental fact of the maximum Tc for M=3 or 4 for a range of 
^\ ■ the parameters considered. That is for moderately small single particle and pair tunneling between successive Cu02 
] layers. As experiments indicate that Tc decreases steadily for M>3 (with the aforementioned exception), we content 
ourselves with calculating Tc for the cases M=l-4, assuming that the trends found also apply for M>4. The complex 
calculations required for M>4 are out of the scope of the present paper. 
We consider the BCS-type Hamiltonian for M CUO2 layers in the unit cell 



i-.k.cr i;k,p 

+ ^^'^ ^k^A-k-T + Y Tk dlj^^^dl_f^ .^di,^k.,idi,k,t ■ 

<il>,k,a <il>,k 

The creation/annihilation operators (i|fe a/^'i;k,a- describe electrons in the z-th Cu02 layer in the unit cell, interacting 
via Vi{k,p), and i =1-M. ei_k = f-i.k~i^ii with the dispersion e^^^ — —2ti{coskx+cosky) — 4t[ cosk^ cosky~2t'-{cos2kx + 
cos2fcy), kx,ky = [— 7r,7r], and the chemical potential fii of the z-th layer. The (in-plane momentum conserving) 
coupling between successive neighbor Cu02 layers < il > for single electron tunneling is ^15i] 

t±,k = to (cos /cj; - COS ky f , (2) 

while for pair tunneling is (proportional to t^ j,, e.g. c.f.[6] and below) 

Tk = Toicoskx - cosfcj,)^ . (3) 
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We diagonalize the single particle kinetic energy part of the Hamiltonian. The original operators di^k,a- are given 
in terms of the new ones Ci^k,tT as dM = Umcm, du — (di,/c,cr, c?2,/c,cr, rfA/.fc.o-) [ill- In. the new basis, we consider 
variational BCS states of the type 

with Ui^kT^i.k the well known coherence factors. 

For M=2 layers ci_2 — {di ± d2)/V2, and the transformed Hamiltonian is 



^^Yl ^''f" ^M-.'y'^^''^''' ^('='P)[cU,t4,-fc,iCi,-p,;Ci,p,t + cl_fc_^c|; _^^^C2,-p4C2^p,t + (1 O 2)] 

+ \ J2 ^fe[cU,tcL-fc,iCi,-fc4Ci,fe,t + 4,fc,tcI,-fc.iC2,-fc4C2,fe,T + (1 ^ 2)] . (5) 

k 

Here ^i.fc = Ek — t±^k, (,2.k = £k + t±,k correspond to the bonding and antibonding states. The two initial layers are 
equivalent, and hence Vi — V2 — V. We also made use of the fact that the only non-zero matrix elements with four 
operators Ci^k,(T (< O >=< 5'|0|5' >) are - c.f. ref. and below for Ui^k,Vi^k and fi^k 

< cIfc,t4-fe4'^J'-P.^'^J'P't >= Vi^kUi,kVj,pUj,p{l - 2/i,fe)(l - 2/j-p) . (6) 

Bearing in mind that uf + vfj, = 1, < u1^,v1^ < 1 , we treat the coherence factors as Ui^k = cos{di^k), 
Vi^k = siii{9i^k)- Minimizing < H >=< ^'|i/|\E' > with respect to 9i^k yields the gap equations 

= -2 Ci^kVi,kU^.k + Ai.kiuij. - vlk] , (7) 

while minimizing < F = H — TS >, with the entropy S = — 2^^ j, [fi^k^T^ fi.k + (1 ^ ./j,/£)lii(l — fi^k)] and the 
temperature T = 1//3, with respect to the thermal factors /i_fc yields 



fi.k -, ftp,. , ; Ei.k — o 2 ' (^) 



The gaps A^^^ are a sum of a diagonal (in the layer index i) part Gi^k and a non-diagonal part gi^k- 

^i,k ^ Gi^k + gi,k , G.i^k ^ - \y^,V{k,p) Vi,pUi,p tanh(/3£;^,p/2) - ^TkVi^kUi^k ta.nh{/3Ei^k/2) , (9) 

p 

9i,k = -lYv{k,p)vj^pUj^p t&nh{pEj^p/2) ~ iTfeW,- fcW,- fc tanh(/3S,,fe/2) , (i, j) = {(1, 2), (2, 1)} . 



2 

p 



Here Aijt = A2.fe. We note that setting i — j and Tk = t±^k — reduces these equations to the usual gap equation 
= -2 ^kVkUk + Ak[ul - vl], with Ak^-Ep V{k,p) VpUp ta.nh{pEp/2) and Ek = VCfc + K ^ ^k/{ul - vl). But, 



note that we do not enforce a relation of the type Ei,k ^ y A^^, 4- S,'^. for M>1. 
The general form of the transformed Hamiltonian is 

H =Y^^,kclk,a(^^,k,a+ ^^ji^^P) Ak,A-k,i'^3-pA(^3:PA ' (10) 

i\k,(T k,p 

The coefficients Wji{k,p) are determined below. Thus we obtain as above the gap equations for Ui_k,Vi.ki fi,k in the 
general form ([T]),®, but now with 

Gi,fe = - ^ Wi,pUi,p tanh(/3£:i^p/2) , g^^fe = ^ -I- i;-,-^pU-,-,p tanh(/3£'j,p/2) . (11) 

P j^'i-.P 



For M=3 layers, S,i^k = £i,fe, C(2,3),fe = (ei.fe + £2,fe T d)/2, with d ^ ^ + ^t\,k^ o. = ei,fc - e2,k- Also 

w,,{k,p) = V^{k,p)[AfMk,p)+j^^Xk,p)] + V2{k,p)AfMk,p) + E,,(k) Skp . (12) 
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The terms ^l™"* correspond to the initial layers m=l-3, with layers 1 and 3 being equivalent: ei^k = £3,k- We have 

A^\k,p) ^ A^^\k,p) = bUk)b',{p) , A^\k,p) = gUk)9]{p) , (13) 
E,j(k) = 2Tfe [ bUk)g^(k) + b^^{k)gUk) ] , 

with bi{k) ^ 1/V2, b2{k) ^ (1 - a/d)ui/ A, bajk) ^ (1 + a/d)u 2/4, giik) = 0, 52(fc) = t±,fcUiM 53(fc) = t^^kU2/d 
and Ml = V2 + (a + rf)V(2ii,fc)', ^2 = ^2 + (« " t^) V(2i±,fc)2. 

For M=4 layers C(i,2).fc = (ei,fe + £2,fe + i±,fc T si)/2, C(3,4),fe = (ei,fe + £2,/c - T S2)/2, with d = ei,fc-e2,fe, si^2 

d2 ± 2 t_L,fc d + 5 J.. Here 

= Vi(A:,p)[4j'(fc,p) + Af{k,p)] + V2{k,p)[A'^\k,p) + Af{k,p)] + E,,{k) 5up . (14) 

The terms ^i™"* correspond to the initial layers 1-4, with layers 1 and 4 and also 2 and 3 being equivalent. Now 

4)(fc,p) = A^{k,p) = bKk)b%p) , A?(fc,p) = A?(A;,p) = <?K^).9,'(P) , (15) 
E,,{k) = 2n [ bi{k)gj{k) + b][k)gUk) + gHk)g^{k) ] , 

with 6i(fc) = ui(d + t_L./c - si)/(4si), 62(fc) = U2{d + <_l./c + si)/(4si), bsik) = U3(-d + ti,fe + S2)/(4s2), biik) = 
M4(c?- ^±,fc + S2)/(4s2), gijk) ^ uit j_^k/{2si ) .g2(fc) = -U2^±,fc/(2si), ff3(fc ) U3<± fc/(252), g4(fc) = M4^±,fc/(2s 2) and 

Ul = y^2+(d + ti,fe+Si)2/(2ti^J, W2 = y^2+(d + ti,fc-Si)2/(2ti^fc), = ^2 + (d - i^.fc + S2)V(24,fe)' "4 - 

^2 + (d~i^,fe-S2)V(2tlfe). 

So far the formalism was quite generic. Specializing to Coulomb repulsion generated positive definite pairing 
potentials [l3] for the cuprates, we consider the realistic non-separable form (and thus harder computationally) 

V,{k,p) = V,{k-p) , VM = Vor sin^ (q^a) sin^ (qya) , Vo^ > , (16) 

which is peaked at (near) Q — (±7r, ±7r) for a = 0.5(0.5 < a < 0.6). This type of potential is well known to generate 



a ci2:2_j,2-wave gap 

mm 



Another relevant issue is the doping imbalance for M>2 layers in the unit cell. Namely, NMR experiments show 



that the outer layers are overdoped with holes compared to the inner layers 17| , in agreement with earlier theoretical 
estimates 0]. We account for this fact by typically considering fj.out = A* ~ 0.06t, /Ui„ = M + 0.06t for both M=3,4, 
with Houtl jJ-in referring to the chemical potentials of the original outer/inner layers and /x the chemical potential for 
the case M=l,2. Assuming a screened electronic interaction, its strength is determined by the susceptibility ^(g, oj) 



12,|l3j. It can be shown that x(g, w = 0) is a decreasing function of ^ for the range of doping considered herein 
12| . Therefore, taking Voi = Vo,out > Vo2 — Vodm reflects also the effect of [lout < Mi" for M=3,4. In our model, 
this potential imbalance is more important than the sheer fj, imbalance. Overall, Vi{k,p) and Tk (c.f. below) mostly 
determine the gaps A^^^ and hence Tc- This simply reflects the fact that both Vi{k,p) and Tk enter on an equal basis 
in the BCS gap equations - c.f. eqs. (9) and (11). The potential has a drastic influence on Tc, as is known from 
standard BCS-EHashberg theory [III [13. 

We give a summary of a few calculations in the 3 figures and the table below, noting that they refer to (ij.2_y2-wave 
solutions of the gap equations above In all cases the energy scale is given hy t = 1, and a — 0.5. Typically values 
0.5 < a < 0.6 yield higher T^'s, without a qualitative change of the results below. 

In the table, an asterisk marks the highest Tc in each case. A maximum Tc for M=3, without resorting to the 
potential imbalance Vo,out > Vo.in {Vi > V2 etc.), can be obtained for a small enough (in magnitude) pair hopping 
To < 0, and this fact is facilitated by smaller values of to - c.f. case A.l and the figures. Otherwise, for Vo,out = Vo.in 
and higher To, M=2 yields the maximum Tc- Tc depends weakly on the single particle hopping to- In ref. [8] it was 
shown that the bare value of to is significantly reduced through the effect of interactions. The sign of to is irrelevant, 
as can be seen from the gap equations. This is not the case for the sign of the pair-coherent term To though (c.f. 
below). To > 0, in general^ does not reproduce the dome- type dependence of Tc on M with a maximum for M=3. 

Given the strong dependence of Tc on Tfc within our model, materials with Y, Ba and Sr intercalants should have 
a much smaller T^ than Ca-intercalated cuprates. Moreover, the in-plane pairing potential is expected to be weaker 
for the former materials for M>1. 

In all, we note that the experimental situation of the Tc increasing for M=l-3 (or 4) and then dropping, is only 
realized for rather small negative values of the pair tunneling amplitude To. E.g. for To = — 0.02t M=2 yields 
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A. t' = -0.35, t" = 0,n = 0.80, n = -0.002 


1. 
2. 
3. 
3. a 
3.b 


to = 0.001 
to = 0.01 
to = 0.001 
to = 0.001 
to = 0.001 


Vol = Vo2 = 4 
Vol = Vo2 = 4 
Vol = Vo2 = 3 

Vol = 3.2, Vo2 = 3 
Vol = 3.5, Vo2 = 3 


Tc{M = 1 - 4) = (5.384, 6.347, *6.424, 6.241)10"^ 
Tc{M = 1 - 4) = (5.384, *6.347, 6.187, 5.852)10"^ 
Tc{M = 1 - 4) = (1.960, *2.943, 2.806, 2.611)10"^ 
T^{M = 3, 4) = (*3.403, 3.240)10^^ 
To{M = 3, 4) = (*4.406, 4.282)10"^ 


B. t' = -0.35, t" = 0, n = 0.80, To = -0.02 


1. 


to = 0.001 


Vol = 4.5, Vo2 = 4 


To{M = 1 - 4) = (5.384, *14.68, 13.02, 9.870)10-^ 


C. t' = -0.3, t" = 0.2, n = 0.85, To = -0.002 


1. 

1. a. 
2. 

2. a. 


to = 0.001 
to = 0.001 
to = 0.01 
to = 0.01 


Vol = Vo2 = 4 

Vol = 4.5, Vo2 = 4 

Vol = Vo2 = 4 

Vol = 4.5, Vo2 = 4 


Tc{M = 1 - 4) = (1.215, *1.700, 1.595, 1.546)10"^ 
To{M = 3, 4) = (*2.264, 2.216)10"^ 

mM = 1 - 4) = (1.215, *1.632, 1.382, 0.9957)10"^ 
rc(M = 3, 4) = (*1.978, 1.627)10"^ 



TABLE I: Tc calculated as a function of the parameters shown. 



the maximum Tc, as opposed to M=3 for To = — 0.002t - c.f. case D.3 below for a maximum for M=4. The 
negative Tk can be understood as follows. Within the frame of second order perturbation theory, we have the 
relation Tu^k — — 2i5_ ^ ''^i^k'^^i^k^i^k^l^k / i^i^k + ^i^k)^ '^i*^ index (0) denoting the initial uncoupled layers, with 

4fe = ^i^Sy+<k- Tu,k < follows from E<^^ + E^^l > 0. 

We give the results of our calculations for five homologuous series of cuprates, which are similar to the cases above. 
We use a restricted set of parameters - 1 = 220 — 245 meV, t' = — 0.35t, t" = 0, to = 0.03t, n = 0.80, which corresponds 
to optimal doping, a = 0.5, To = — 0.0044t — (— 0.0018)t, and Voi = 2.7t — 4.2t - to demonstrate the fitting capacity of 
our model. We obtain exact matches with the experimental Tc values, given in degrees K and taken from ref. We 
also give the ratios R=max{Ai^fe(T = 0)}/Tc, which turn out to be in the range 0.71-1. These values are too small 
compared to the actual experimental ones, pointing to the limitations of the BCS description for the cuprates. Along 
the same line, the Voi values are in the intermediate coupling regime, i.e. beyond the strict limit of applicability of the 
BCS framework. Note that the opposite potential imbalance Vi < V2 would yield the same Tc's with somewhat higher 
values of V2 for M=3,4 than the Vi values below. The values of t and Vol are adjusted so as to yield Tc(M=l). Then, 
To is chosen so as to yield Tc(M=2). Subsequently, Vo2(M = 3,4) = Voi{M = 1) = Voi(M = 2), and Ki(M = 3,4) 
are chosen so as to yield Tc(M=3,4). 

D.l Bi2Sr2CaM-iCuM02M+4+5 with t = 245 meV and To = -0.0044t, Ki = V02 = 2.7t. For M=3 we take 
Vol = 3.5t, thus obtaining re(M=l-3)=36, 90, 110 K and R(M==l-3)==0.769, 0.713, 0.867. 

D.2 TlBa2CaM-iCuM02M+3+5 with t = 229 meV and To = -0.0044t, Voi = V02 = 3t. For M=3 we take 
Vol = 3.85t and for M=4 Voi = 4.16t, thus obtaining rc(M=l-4)=52, 107, 133, 127 K and R(M=l-4)=0.765, 0.781, 
0.973, 1.034. 

D.3 Tl2Ba2CaAf-iCuA/02Af+3+5 with t = 220 meV and To = -0.0024t, Voi = V02 = 3t. For M=3 we take 
Vol = 3.67f and for M=4 Voi = 4.2t, thus obtaining rc(M=l-4)=50, 80, 110, 122 K and R(M=l-4)=0.765, 0.755, 
0.909, 1.039. This series has maximum Tc for M=4 - c.f. 0. 

D.4 Tl2Ba2CaA/_iCuA/02M+4+5 with t = 221 meV and To = -0.0018t, Voi = V02 = 3.5t. For M=3 we take 
Vol = 3.76* and for M==4 Voi = 3.83t, thus obtaining Tc(M=l-4)==90, 115, 125, 116 K and R(M=l-4)=0.795, 0.845, 
0.934, 0.976. 

D.5 HgBa2CaA,/_iCuAf02A/+2+A- with t = 238.5 meV and To ^ -0.0021t, Ki = K2 = 3.5t. For M=3 we take 
Vol = 3.75t and for M==4 Voi = 3.89t, thus obtaining Tc{M^\-A)^97 , 127, 135, 129 K and R(M=l-4)=0.765, 0.851, 
0.932, 0.979. 

In summary, higher Tc in multilayer copper oxides arises from a relatively increased (in magnitude) pair interlayer 
coupling Tfe < 0, combined with a substantial strength of the repulsive Coulomb interaction between the electrons. 
It seems that all these ingredients are present in the Ca-intercalated materials. This picture should possibly be 
complemented by the Cu02 lattice symmetry effects [31 . The BCS approach, though of limited applicability in the 
cuprates, offers relevant insight. We give an interpretation of the negative Tk in the BCS framework. 

Stamatis Stamatiadis provided essential programming assistance. The author enjoyed discussions with Josheph 
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FIG. 1: Tc (in units of t) as a function of interplane hopping for n — 0.85, t' ■- 
layer doping imbalance is taken into account for M=3,4. Squares : M=2, x's 
for To — 0. All lines are guides to the eye. 
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FIG. 2: Same as in fig. 1. Tc versus to for To = -0.002 
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FIG. 3: Same as in fig. 1. Tc versus (-To) for to = 10 . Notice the log scale for Tc. 



